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We show how the squeezing of light can lead to the formation of topological states. Such states 
are characterized by non-trivial Chern numbers, and exhibit protected edge modes which give rise 
to chiral elastic and inelastic photon transport. These topological bosonic states are not equivalent 
to their fermionic (topological superconductor) counterparts and cannot be mapped by a local 
transformation onto topological states found in particle-conserving models. They thus represent a 
new type of topological system. We study this physics in detail in the case of a Kagome lattice model, 
and discuss possible realizations using nonlinear photonic crystals or superconducting circuits. 


I. INTRODUCTION 

There has been enormous interest in trying to replicate 
the physics of topological electronic phases in a variety 
of bosonic systems, including cold atoms [1], photonic 
systems [2] and more recently phononic systems nu. 
Photonic analogues include quantum-Hall like states in¬ 
duced through the introduction of synthetic gauge fields 
mm, phases which are analogous to time-reversal in¬ 
variant topological insulators [ME], Floquet topolog¬ 
ical insulators [201 ED and even Major ana-like modes 
[22] . Experimental studies of such phases have made sig¬ 
nificant progress [19], 21, f23H26] . In addition to being 
of fundamental interest, these topological photonic and 
phononic phases could have practical utility, as they pro¬ 
vide disorder-protected edge modes that could be used 
for chiral light and sound propagation. 

Despite this intense activity, most works on topological 
bosonic states amount in the end to replicating a well- 
known fermionic single-particle Hamiltonian with bosons; 
as the topological properties are a function of the result¬ 
ing single-particle eigenstates, particle statistics play no 
crucial role, except perhaps in the methods used for prob¬ 
ing the system. As we now discuss, this simple correspon¬ 
dence will fail if the particle number is not conserved. 

Consider the most general quadratic Hamiltonian de¬ 
scribing bosons on a lattice which respects the discrete 
translational invariance of the lattice, but which does not 
conserve particle number: 

H = X e «[ k ] 6 tA > ™+ V (v«'[k]&b/-k,n' + h - c ) 

k,n k ,n,n' 

(i) 

The first term describes a non-interacting band- 
structure, where k runs over momenta in the first Bril- 
louin zone, and n labels the bands. The remaining terms 
correspond to parametric driving or two-mode squeezing 
terms. As we discuss below, they can be controllably 
realized in a number of different settings, with the 5k, m 
operators describing "real" particles (i.e. photons in a 
cavity lattice), and not quasiparticles defined above some 



Figure 1. Setup figure: (a) An array of nonlinear cavities 
forming a kagome lattice, (b) Photons hop between nearest- 
neighbor sites with rate J. Each cavity is driven parametri¬ 
cally leading to the creation of photon pairs on the same lat¬ 
tice site (rate v on ) and on nearest-neighbor sites (rate v 0 r)- A 
spatial pattern of the driving phase is imprinted on the para¬ 
metric interactions, breaking the time-reversal symmetry (but 
preserving the C 3 rotational symmetry). 

effective condensate. While superficially similar to pair¬ 
ing terms in a superconductor, these two-mode squeezing 
terms have a profoundly different effect in a bosonic sys¬ 
tem, as there is no limit to the occupancy of a particular 
single-particle state. They can give rise to highly entan¬ 
gled ground states, and even to instabilities. 

Given these differences, it is natural to ask how anoma¬ 
lous "pairing" terms can directly lead to topological 
bosonic phases. In this work, we study the topological 
properties of 2D systems described by Eq. Q , in the case 
where the underlying particle-conserving band structure 
has no topological structure, and where the parametric 
driving terms do not make the system unstable. We show 
that the introduction of particle non-conserving terms 
can break time-reversal symmetry (TRS) in a manner 
that is distinct from having introduced a synthetic gauge 
field, and can lead to the formation of bands having a 
non-trivial pattern of (suitably defined) quantized Chern 
numbers. This in turn leads to the formation of protected 
chiral edge modes: unlike the particle conserving case, 
these modes can mediate a protected inelastic scattering 
mechanism along the edge (i.e. a probe field injected into 
the edge of the sample will travel along the edge, but 
emerge at a different frequency). In general, the topo- 
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logical phases we find here are distinct both from those 
obtained in the particle-conserving case, and from those 
found in topological superconductors. We also discuss 
possible realizations of this model using a nonlinear pho¬ 
tonic crystal or superconducting microwave circuits. 

Note that recent works have explored topological fea¬ 
tures of bosonic quasiparticles in condensed phases; in¬ 
teractions with the condensate treated at the mean-field 
level yield a Hamiltonian having the general form of 
Eq. 0. These include a study of a magnonic crystal m, 
as well as general Bose-Einstein condensates in ID [28] 
and in 2D [29] . Unlike our work, these studies did not ex¬ 
plore the general role that tunable squeezing terms play 
in yielding topological states. Further, in our case Eq. (1) 
describes the M real M particles of our system, not quasipar¬ 
ticles defined above some background; this makes detec¬ 
tion and potential applications much easier. Most impor¬ 
tantly, the topologically-protected inelastic scattering we 
describe is absent in those different settings. 

II. RESULTS 
A. Kagome lattice model 

For concreteness, we start with a system of bosons on 
a kagome lattice (see Fig0i 

Ho = yy o ajaj - J ^2 «j«r (2) 

j <jJ') 


be achieved by a corresponding variation of the driving 
phase of the pump. Note that we are working in a rotat¬ 
ing frame where this interaction is time-independent, and 
thus u)q should be interpreted as the detuning between 
the parametric driving and the true on-site (cavity) fre¬ 
quency Ci; cav (i.e. cjo = u; cav — cjl/ 2, where the parametric 
driving is at a frequency ujl). The parametric driving can 
cause the system to become unstable; we will thus require 
that the on-site energy uj 0 is always sufficiently large that 
each parametric driving term is sufficiently non-resonant 
to ensure stability. 

For a generic choice of phases in the parametric driv¬ 
ing Hamiltonian of Eq. (J3|, it is no longer possible to 
find a gauge where H = Hq + Hl is purely real when 
expressed in terms of real-space annihilation operators: 
hence, even though the hopping Hamiltonian Hq corre¬ 
sponds to strictly zero flux, the parametric driving can 
itself break TRS. In what follows, we will focus for sim¬ 
plicity on situations where time-reversal and particle- 
conservation are the only symmetries broken by the para¬ 
metric driving: they will maintain the inversion and C3 
rotational symmetry of the kagome lattice. We will also 
make a global gauge transformation so that v 0 r is purely 
real, while v on = In this case, the only possible 

choices for the <f phases have the form (4>A,f>B,4>c) = 
{(^abAbcAca) = ±(0,5,25) with S = 27rra*,/3, m v is 
an integer and is the vorticity of the parametric pumps. 

B. Gap opening and non-trivial topology 


(we set h — 1). Here, we denote by dj the photon anni¬ 
hilation operator associated with lattice site j, where the 
vector site index has the form j = (ji,j2js)- ji,j 2 £ Z 
labels a particular unit cell of the lattice, while the index 
s = A,B,C labels the element of the sublattice, (j, j') 
indicates the sum over nearest neighbors, and J is the 
(real valued) nearest neighbor hopping rate; c^o plays the 
role of an on-site energy. As there are no phases associ¬ 
ated with the hopping terms, this Hamiltonian is time- 
reversal symmetric and topologically trivial. We chose 
the kagome lattice because it is directly realizable both 
in quantum optomechanics [5j and in arrays of super¬ 
conducting cavity arrays puma; it is also the simplest 
model where purely local parametric driving can result 
in a topological phase. 

We next introduce quadratic squeezing terms to this 
Hamiltonian that preserve the translational symmetry of 
the lattice and that are no more non-local than our orig¬ 
inal, nearest-neighbor hopping Hamiltonian: 


Hr= -- 




at a | 


-v 0 fi 2^ 1 
<j»j') 


'44 


( 3 ) 

Such terms generically arise from having a nonlinear in¬ 
teraction with a driven auxiliary pump mode on each 
site (which can be treated classically). As we discuss be¬ 
low, the variation in phases in Hl from site to site could 


Hq is the standard tight-binding kagome Hamilto¬ 
nian for zero magnetic field, and does not have band 
gaps: the upper and middle bands touch at the sym¬ 
metry point T = (0,0), whereas the middle and lower 
bands touch at the symmetry points K = (27r/3, 0) and 
K' = (7r/3, 7 t/\/ 3) where they form Dirac cones [see 
Fig- Ha)]. 

Turning on the pairing terms, the Hamiltonian H = 
Hq + Hl can be diagonalized in the standard manner as 
H = En[k\Pl k fJ rhk - where the /3„ jk are canonical 
bosonic annihilation operators determined by a Bogoli- 
ubov transformation of the form (see Appendix [A] for 
details): 

/^n,k == ^n,k[s]flk )S — ^n,k[s]®—k,s- (4) 

s=A,B,C 

Here, dk, s are the annihilation operators in quasi¬ 
momentum space, and n = 1,2,3 is a band index; 
we count the bands by increasing energy. The pho¬ 
tonic single-particle spectral function now shows reso¬ 
nances at both positive and negative frequencies, ±E n [k \, 
corresponding to “particle”- and “hole”-type bands, see 
Fig. [2|d). Because of the TRS breaking induced by the 
squeezing terms, the band structure described by E n [ k] 
now exhibits gaps, see Fig. §b); further, for a finite 
sized system, one also finds edge modes in the gap, see 
Fig. Hp). 
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Figure 2. Band structure: (a-b) 3D plots of the bulk band structure for J — 0.02o;o- The hexagonal Brillouin zone is also 
shown, (a) In the absence of parametric driving, neighboring bands touch at the rotational symmetry points K and K' and T. 
(b) The parametric driving opens a gap between subsequent bands. For the chosen parameter (y on = —0.085 and v 0 r = 0.22) 
there is a global band gap between the second and third band, (d) Hole and particle bands, ±Em[k x \, in a strip geometry 
[sketched in (c)]. The line intensity is proportional to the weight of the corresponding resonance in the photon spectral function, 
see Appendix | e| The edge states localized on the right (left) edge, plotted in green (blue), have positive (negative) velocity. 


The above behaviour suggests that the parametric 
terms have induced a non-trivial topological structure in 
the wavefunctions of the band eigenstates. To quantify 
this, we first need to properly identify the Berry phase 
associated with a bosonic band eigenstate in the pres¬ 
ence of particle non-conserving terms. For each k, the 
Bloch Hamiltonian H k corresponds to the Hamiltonian of 
a multi-mode parametric amplifier. Unlike the particle- 
conserving case, the ground state of such a Hamiltonian is 
a multi-mode squeezed state with non-zero photon num¬ 
ber; it can thus have a non-trivial Berry’s phase asso¬ 
ciated with it when k is varied, see Appendix [B] The 
Berry phase of interest for us will be the difference of 
this ground state Berry phase and that associated with 
a single quasiparticle excitation. One finds that the re¬ 
sulting Berry connection takes the form 

An = i(k,n|^V fc |k,n), (5) 

Here, the 6-vector of Bogoliubov coefficients |k, n) = 

(^n,k \A\ , [-5] 5 ^n,k [ 6 ] •> ^n,k [A\ , V n \ c [-B], U n? k [ 6 ]) plays 

the role of a singe-particle wavefunction, and a z acts 
in the ’particle-hole-space’, associating +1 to the u - 
components and —1 to the u-components, see Appendix 
|B| for details. These effective wavefunctions obey the 
symplectic normalization condition 

(k,n|<j z |k', n') = ^< )k H«n',k'[s] - v* >k [s]tv,k'[s] 

= ^k,k (6) 

Having identified the appropriate Berry connection for 
a band eigenstate, the Chern number for a band n is then 
defined in the usual manner: 

C n = T [ (V X A n ) ■ Z (7) 

Jbz 

This definition agrees with that presented in Ref. m 
and (in ID) Ref. J28|; standard arguments |27j show that 


the C n are integers with the usual properties. We note 
that, as for superconductors, breaking the U( 1) (particle- 
conservation) symmetry remains compatible with a first- 
quantized picture after doubling the number of bands. 
The additional ’hole’ bands are connected to the stan¬ 
dard ’particle’ bands by a particle-hole symmetry; see 
Appendix [A[ In bosonic systems, the requirement of sta¬ 
bility implies that particle and hole bands can not touch. 
Thus, the sum of the Chern numbers over the particle 
bands (with E > 0) must be zero, and there cannot be 
any edge states with energies below the lowest particle 
bulk band (or in particular, at zero energy); see Appendix 

IB 

In the special case where we only have onsite paramet¬ 
ric driving (i.e. z/ Q ff = 0, v on ^ 0), the Chern numbers 
can be calculated analytically (see Appendix |C|). They 
are uniquely fixed by the pump vorticity. If m u = 0, 
we have TRS and the band structure is gapless, while 
for m v = ±1, 6 = (=Fl,0, ±1). This set of topologi¬ 
cal phases also occurs in a particle-number conserving 
model on the kagome lattice with a staggered magnetic 
field, i.e. the Oghushi-Murakami-Nagaosa (OMN) model 
of the anomalous quantum Hall effect [301 EQ. 

In the general case, including offsite parametric driv¬ 
ing, entirely new phases appear. We have computed the 
Chern numbers for that case numerically, using the ap¬ 
proach of Ref. [32] • In Fig. |3|a), we show the topo¬ 
logical phase diagram of our system, where J/ujq and 
rrijy are held fixed, while the parametric drive strengths 
^on, ^off are varied. Different colors correspond to differ¬ 
ent triplets C = ( 61 , 62 , 63 ) of the band Chern num¬ 
bers, with gray and dark-gray corresponding to the two 
phases already present in the OMN model. Strikingly, 
a finite off-diagonal coupling v 0 r generates a large vari¬ 
ety of phases which are not present in the OMN model. 
The border between different topological phases repre¬ 
sent topological phase transitions, and correspond to pa¬ 
rameter values where a pair of bands touch at a particular 
symmetry point; we discuss this further below. 
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Figure 3. (a) Topological phase diagram for the parametrically 
driven kagome lattice model. The y (tc) axis corresponds to 
the strength of the onsite parametric drive v on (off-site para¬ 
metric drive z'off), and different colours correspond to dif¬ 
ferent triplets C = (Ci, C 2 , C 3 ) of Chern numbers for the 
three bands of the model. Note that only the gray and dark- 
gray phases are found in the particle-conserving version of our 
model with a staggered field. We have fixed the hopping rate 
J/ujo — 0.02, and the vorticity of the pump mn u — 1. (b) 
Same phase diagram, but now plotted in terms of the effec¬ 
tive flux <f> and effective parametric drive v experienced by a 
quasiparticles. See text for details. 


C. Effective model 


To gain further insight into the structure of the topo¬ 
logical phases found above, it is useful to work in a 
“dressed-state" basis that eliminates the local parametric 
driving terms from our Hamiltonian. We thus first diago¬ 
nalize the purely local terms in the Hamiltonian; for each 
lattice site j we have 


H\ = cd 0 a a. - - 


1 


2 L 


z/ on e^ j djdj + h.c. 


= &a\a .. ( 8 ) 


Here uj = and the annihilation operators 

dj are given by a local Bogoliubov (squeezing) transfor¬ 
mation dj = e^e - ^/ 2 ^coshraj — e l ^e L(pv sinh raj^, 
where the squeezing factor r is 


r = 



^0 + ^on 
CJq lS 0 n 


(9) 


On a physical level, the local parametric driving terms 
attempt to drive each site into a squeezed vacuum state 
with squeeze parameter r; the dj quasiparticles corre¬ 
spond to excitations above this reference state. Note that 
we have included an overall phase factor in the definition 
of the dj which will simplify the final form of the full 
Hamiltonian. 

In this new basis of local quasiparticles, our full Hamil¬ 
tonian takes the form 


j (j>i) 


+h ' c - 


<j,i> 


( 10 ) 


The transformation has mixed the hopping terms with 
the non-local parametric terms: The effective counter¬ 
clockwise hopping matrix element is 


Jji = Je l 


a 3iS/2 


2 J cos sinh 2 r + z/ Q ff sinh 2r cos ( ^ 




(ii) 


and the magnitude of the effective non-local parametric 
driving is 

\v\ = \iy 0 ^e~ l ^ s ^ 2+cp ^ + 2 z/ 0 ff cos (£/2 + (p u ) sinh 2 r 

+ J sinh2r cos(5/2)|. ( 12 ) 

Note that the phase of v can be eliminated by a global 
gauge transformation, and hence it plays no role; we thus 
take v to be real in what follows. 

Our model takes on a much simpler form in the new 
basis: the onsite parametric driving is gone, and the non¬ 
local parametric driving is real. Most crucially, the ef¬ 
fective hoppings now can have spatially-varying phases, 
which depend both on the vorticity of the parametric 
driving in Hl (through S), and the magnitude of the on¬ 
site squeezing (through r). In this transformed basis, the 
effective hopping phases are the only route to breaking 
TRS. If the parametric terms were not present, the com¬ 
plex phases would correspond in the usual manner to a 
synthetic gauge field (i.e. the effective flux piercing a 
triangular plaquette would be = 3arg J). Our model 
has thus been mapped on to the standard OMN model 
for the anomalous quantum Hall effect, with an addi¬ 
tional (purely real) nearest-neighbour two-mode squeez¬ 
ing interaction. However, we note that strictly speaking 
4> can not be interpreted as a flux in the presence of the 
additional parametric terms: a ’flux’ of 2tt can not be 
anymore eliminated by a gauge transformation because 
the complex phases reappear in the parametric terms. In 
that case, only a periodicity of 67 r in <f> is retained, since 
that corresponds to having trivial hopping phases of 27 t. 

Understanding the topological structure of this trans¬ 
formed Hamiltonian is completely sufficient for our pur¬ 
poses: one can easily show that the Chern number of a 
band is invariant under any local Bogoliubov transforma¬ 
tion, hence the Chern numbers obtained from the trans¬ 
formed Hamiltonian in Eq. © will coincide exactly with 
those obtained from the original Hamiltonian in Eq. ©>• 
We thus see that the topological structure of our system 
is controlled completely through only three dimension¬ 
less parameters: the ’flux’ <f> associated with the hopping 
phase, the ratio \v/J\ and the ratio uj/\J\. 

The topological phase diagram for the effective model 
is shown in Fig. [3](b) . Again, one sees that as soon as 
the effective non-local parametric drive v is non-zero, 
topological phases distinct from the standard (particle- 
conserving) OMN model are possible. The sign of the 
parametric pump vorticity m v determines the sign of the 
effective flux 4>, c.f. Eq. |TT] As such, the right half of 
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Fig. [3jb) (corresponding to $ > 0) is a deformed version 
of the phase diagram of the original model for pump vor- 
ticity m v = 1, as plotted in Fig. |3ja). Changing the sign 
of m v (and hence 4>) simply flips the sign of all Chern 
numbers, see Appendix 0 

Our effective model provides a more direct means for 
understanding the boundaries between different topolog¬ 
ical phases. Most of these are associated with the cross¬ 
ing of bands at one or more high-symmetry points in the 
Brillouin zone; this allows an analytic calculation of the 
phase boundary (see AppendixjC]). Perhaps most striking 
in Fig. |3|b) is the horizontal boundary (labelled A4), oc¬ 
curring at a finite value of the effective offsite parametric 

drive, v « yf~JCb. This boundary is set by the closing of a 
band gap at the M points; as these points are associated 
with the decoupling of one sublattice from the other two, 
this boundary is insensitive to the flux 4>. Similarly, the 
vertical line labeled T denotes a line where the system 
has TRS, and all bands cross at the symmetry points K , 
K' and T. The case of zero pump vorticity m u = 0 (not 
shown) is also interesting. Here, the effective flux <£> de¬ 
pends on the strength of the parametric drivings, but is 
always constrained to be 0 or 37r. This implies that the 
effective Hamiltonian has TRS, even though the original 
Hamiltonian may not (i.e. if Im v q r ^ 0, the original 
Hamiltonian does not have TRS). For m v = 0, the para¬ 
metric drivings do not open any band gap and the Chern 
numbers are not well defined. 


D. Edge states and transport 

Despite their modified definition, the Chern numbers 
associated with our Bogoliubov bands still guarantee the 
existence of protected chiral edge modes in a system 
with boundaries via a standard bulk-boundary corre¬ 
spondence, see Appendix [Dj These states can be used 
to transport photons, by exciting them with an auxiliary 
probe laser beam which is focused on an edge site and 
at the correct frequency. The lack of particle-number 
conservation manifests itself directly in the properties of 
the edge states: along with the standard elastic trans¬ 
mission they can also mediate inelastic scattering pro¬ 
cesses. In terms of the original lab frame, light injected 
at a frequency uo p can emerge on the edge at frequency 
ml — Up where uol is the frequency of the laser paramet¬ 
rically driving the system. This is analogous to the idler 
output of a parametric amplifier. Here, both signal and 
idler have a topologically protected chirality. 

Shown in Fig. [4] are the results of a linear response 
calculation describing such an experiment, applied to a 
finite system with corners. We incorporate a finite pho¬ 
ton decay rate k in the standard input-output formalism, 
see Appendix |E] for details. Narrow-band probe light in¬ 
side a topological band gap is applied to a site on the 
edge, and the resulting inelastic transmission probabili¬ 
ties to each site on the lattice are plotted, see Fig. |4ja). 
One clearly sees that the probe light is transmitted in a 



momentum k at <£ 

Figure 4. Topologically protected transport in a finite sys¬ 
tem: (a) A probe beam at frequency u p inside the bulk band 
gap is focused on the edge of a finite sample. The proba¬ 
bility map of the light transmitted inelastically at frequency 
ujl — oj p clearly shows that the transport is chiral, (b) The 
elastic and inelastic transmission probability to a pair of sites 
along the edges [indicated in red in (a)] is plotted in blue and 
green, respectively. A cut through the bulk bands is shown 
to the left. (c,d) Sketch of the relevant scattering processes 
and energy scales. The inelastic (elastic) transmission has a 
larger rate when the light is injected in the hole ( particle) band 
gap. Parameters: J/w o = 0.02, u 0 n/wo = 0.4, u 0 e/ujo = 0.02, 
k/cjq — 0.001. In (a) uj p — ujl/^ — 0.95o;o 


uni-directional way along the edge of the sample, and is 
even able to turn the corner without significant backscat- 
ter. The corresponding elastic transmission [not shown] 
is also chiral and shows the same spatial dependence. 
In Fig. gb) we show the elastic and inelastic transmit- 
tions to the sites indicated in red (rescaled by the overall 
transmission, 1 — R where R is the reflection probability 
at the injection site) as a function of the probe frequency 
uj p . By scanning the laser probe frequency one can sepa¬ 
rately address particle and hole band gaps. The relative 
intensity of the inelastic scattering component is highly 
enhanced when the probe beam is inside a hole band gap, 
see also the sketches in Fig.^c-d). When the parametric 
interaction between the a quasiparticles is negligible, the 
ratio of elastic and inelastic transmissions depends only 
on the squeezing factor r, [c. f. Eq. ([9|], see Appendix [E] 


E. Physical realization 

Systems of this type could be implemented in 2D pho¬ 
tonic crystal coupled cavity arrays fabricated from non¬ 
linear optical materials [33H35] . The array of optical 
modes participating in the transport would be supple¬ 
mented by pump modes (resonant with the pump laser 
at twice the frequency). One type of pump mode could be 
engineered to be spatially co-localized with the transport 
modes (z/ on processes), while others could be located in- 
between {y off). The required periodic phase pattern of 
the pump laser can be implemented using spatial light 
modulators or a suitable superposition of several laser 
beams impinging on the plane of the crystal. Optome- 
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chanical systems offer another route towards generating 
optical squeezing terms PH EH, via the mechanically in¬ 
duced Kerr interaction, and this could be exploited to 
create an optomechanical array with a photon Hamilto¬ 
nian of the type discussed here. Alternatively, these sys¬ 
tems can be driven by two laser beams to create phononic 
squeezing terms [38]. A fourth alternative consists in su¬ 
perconducting microwave circuits of coupled resonators, 
where Josephson junctions can be embedded to introduce 
and higher-order nonlinearities, as demonstrated in 

mm- 


III. CONCLUSIONS AND OUTLOOK 


In this work, we have shown how tunable squeezing 
interactions in a photonic system are flexible tools that 
allow the creation of new kinds of topological bosonic 
phases. We further demonstrated that the protected edge 
channels supported by such phases allow uni-directional 
elastic and inelastic coherent photon transport. Our work 
opens the door to a number of interesting new directions. 
On the more practical side, one could attempt to ex¬ 
ploit the unique edge states in our system to facilitate 
directional, quantum-limited amplification. On the more 
fundamental level, one could use insights from the corre¬ 
sponding disorder problem m and attempt to develop 
a full characterization of particle non-conserving bosonic 
topological states that are described by quadratic Hamil¬ 
tonians. This would then be a counterpart to the classi¬ 
fication already developed for fermionic systems [42] . 

VP, CB, and FM acknowledge support by an ERC 
Starting Grant OPTOMECH, by the DARPA project 
ORCHID, and by the European Marie-Curie ITN net¬ 
work cQOM. MH and AC acknowledge support from 
NSERC. 


Appendix A: Details of the calculation of the band 
structure 


1. Bogoliubov transformation and first-quantized 
picture 


Here, we show how to find the Bogoliubov transfor¬ 
mation Eq. Q which diagonalizes the Hamiltonian H 
defined in Eqs. It is convenient to switch to a 

first quantized picture by casting the second quantized 
Hamiltonian H in the form 


Me (all '•« 




where dk = (dkA, dkc) T and d(k) is the Bogoliubov 
de Gennes Hamiltonian. By plugging the Bogoliubov 
ansatz Eq. @ into the Heisenberg equation of motion 


for the normal modes ladder operators k , 

one immediately finds the generalized eigenvalue problem 

ft(k)|k n ) = E n [k\a z \k n ). (Al) 

Here, the matrix a z is (minus) the identity for the (hole) 
particle sector, 


The bosonic quadratic Hamiltonian H describes the dy¬ 
namics linearized around a stable classical solution if it 
is possible to find its normal mode decomposition. In 
other words, for all k the eigenvalue problem Eq. © 
should have a set of three orthonormal (particle-like) so¬ 
lutions, c. f. Eq. §• A sufficient and necessary criterion 
for stability is that all eigenvalues of & z h(k) are real, see 
discussion below. 


2. Particle-hole symmetry 


The symplectic eigenvalue problem © has an embed¬ 
ded particle-hole symmetry: for any particle-like solution 
with momentum k, energy E n [ k], and wavefunction |k n ) 
there is a hole-like solution with momentum — k, energy 
—E n [ k], and wavefunction 

| - k_ n ) = a x JC\k n ) = C\k n ). 

(notice that the hole-like solutions have negative length 
(k_ n \a z \k_ n ) = —1). Here, operator /C gives the com¬ 
plex conjugate of the wavefunction and a x exchanges par¬ 
ticles and holes, 


O X 




In other words, the Bogoliubov de Gennes Hamiltonian 
has the generalized symmetry C^h(k)C = — h(— k) where 
the charge conjugation operator C is anti-unitary and 
C 2 = Hq. Thus, our system represents the Bosonic ana¬ 
logue of a superconductor in the Class D of the standard 
topological classification. 


3. Numerical calculation of the band structure 


From Eqs. ( 2|3 ) one can immediately derive the explicit 
expression of the Bogoliubov de Gennes Hamiltonian 


H k) = 


w 0 - Jr(k) 

, d( k ) 


M k) \ 

uj 0 - Jf(k)J 


(A2) 


where Hl = —v on K mi ' — i ,, oS (A mi 'T + fA 1 " 1 ') (m^ is 
the pump vorticity). The matrices f(k) and A are 
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single-particle operators acting on the sub-lattice de¬ 
grees of freedom. They have matrix elements: A a a = 
1, A bb = e i2w / 3 , A CC = e~ 2i ^ 3 , t ba = t* ab = 
1 + e lk ai , t CB = t* bc = 1 + e* k ' aa , t ac = t* ca = 
1 + e* k “ 3 where ai = (—1,— \/3), = (2,0), and 

03 = (-1,V3) are lattice vectors. All other matrix ele¬ 
ments are zero. Notice that A raises the quasi-angular 
momentum by one unit: A|p, m) = \p, m + 1) where 
|p, m) = (1, exp[/2ra7r/3], exp[—i2ra7r/3], 0,0, 0) and like¬ 
wise for the holes. Here, m = 0 is the vortex free eigen¬ 
state, m = 1 has a vortex, and m = — 1 an anti-vortex 
(m is defined modulo 2). Thus, a hole with quasimomen¬ 
tum m is converted into a particle with quasi-angular 
momentum m + m v . In other words, a pair of down- 
converted photons have quasi-angular momenta — m and 
m^+m, respectively. The additional quasimomentum m u 
is provided by the pump photons. In order to write the 
off-diagonal parametric interaction compactly in terms 
of the quasi-angular momentum raising operator A we 
have introduced the rescaled off-diagonal parametric cou¬ 
pling z/ ff = e _2(5 / 2 ^off for pump vorticity m u = ±1 and 
z/ ff = ^off/2 for pump vorticity m v = 0. 

For the effective model Eq. ®> the Bogoliubov de 
Gennes Hamiltonian reads 

Mk) = (“ “ M <*;*> . -^< k) J . (A3) 

\ — vt ( k) uj — Jr(k, —4>)/ 


nian 

Hl{ k) = - on + v'oA T m + T m - m J^rr^llc mv _ m 

m 

+H.c., for k = K,K', 

hl {o = -{ y] Kn + v' 0 ff (r m + 

m^m u 

— {Von/2 + H.C. 


where ak, m is the creation operator of an excitation with 
quasimomentum k (k = T, K, K r ) and quasi-angular mo¬ 
mentum m. Moreover, r m indicates the correspond¬ 
ing eigenvalue of the hopping matrix f: ro(T) = 4, 
n(T) = r_!(r) = n(K) = = -2, TO(10 = 

Ti(AT') = to(K') = r-i(K) = 1. By diagonalizing the 
squeezing Hamiltonian we find the general expression for 
the eigenvalues 


77 , T I'm 1~m—m u . 

J ~ r 


UJQ ~ J 


I'm T 7 ~m—m v 


|^on T ^ 0 ff ("Tn H - 7 ~m—r 


-I 1/2 


(A4) 


From an analogous calculation we obtain the spectrum 
of the effective model at the rotation symmetry points 


En 



2 ~2 
T V 


1/2 


(A5) 


Here, we have introduced the hopping matrix f (k, 4>) 
in the presence of the synthetic magnetic field flux 4>. 
It has the following non zero matrix elements: tba — 
t*ab = e i$/3 (l + e ikttl ), f CB = t* BC = e^ 3 (l + e ik “ 2 ), 
t A c — ^ca = e* $ / 3 (l + e ,k ' a:i ). In the most general case, 
we calculate the band structure and the single-particle 
wavefunctions by diagonalizing the 6x6 matrix cr z h{ k) 
numerically. 


4. Analytical calculation of the band structure 
close to the symmetry points 

One can gain much insight on the array dynamics, in¬ 
cluding the stability requirements and the array topol¬ 
ogy, by calculating analytically the band strucure at the 
symmetry points. This is a particularly easy task at the 
rotational symmetry points T, A, and K '. There, the 
hopping matrix f is diagonal in the basis of the quasi- 
angular momentum eigenstates. Thus, in this basis the 
Hamiltonian becomes block diagonal with 2-dimensional 
blocks. Each block is described by a two-mode squeezing 
Hamiltonian, except for the quasi-momentum k = T and 
m = —m Ul when it is a single mode squeezing Hamilto¬ 


where f m (k, <f>) are the eigenvalues of f: fo(T) = 
4cos[<f>/3], r± i(T) = 4cos[27r/3 =p <E>/3], ri(K) = 
T-xiK') = —2cos[4>/3], n(K') = t 0 (K) = 2cos[tt/3 - 
4>/3], T-i(K) = ro(K') = 2cos[tt/3 + 4>/3]. 

5. Stability analysis 

The system is stable when all eigenvalues of <r z h( k) 
are real. If all eigenenergies of the unperturbed Hamil¬ 
tonian Hq have the same sign the parametric interaction 
is off-resonant and the system is stable if the parametric 
couplings are below a finite threshold. On the contrary, 
if the unperturbed band touches the zero-energy axis, 
the parametric interaction is resonant for the zero energy 
modes leading to an instability for any arbitrarily small 
value of the coupling. Thus, the parametric instability 
sets an upper limit to the hopping J. For concreteness, 
we consider a positive onsite energy ujo (corresponding 
to a red detuned drive). In this case, all eigenenergies 
of Ho are positive if J < cjo/ 4. In this case, the system 
is stable for sufficiently small values of the parametric 
couplings z/ on and v q r. Nevertheless, the threshold of an 
instability is reached as soon as the lowest eigenenergy of 
a particle-type band becomes zero for a finite value of the 
parametric couplings v on and v Q &. For the parameters of 
Fig-! the lowest band touches the zero energy axis at 
the T point. Thus, we can find an analytical expression 








for the instability threshold usi ng the so lutions at the ro¬ 
tational symmetry points Eqs. ( A4|A5 ). In Fig. 3(a) the 
state with zero energy at the border of the unstable re¬ 
gion has zero quasimomentum and an anti-vortex. Thus, 
the instability threshold is given by setting A_i(T) = 0. 
In Fig. 3(b) the state with zero energy has zero quasi¬ 
momentum and quasi-angular momentum. By setting 
E 0 (T) = 0, we find a simple expression for the instability 
threshold, v = uj /4 — Jcos[4>/3]. 


Appendix B: Details of the definition and properties 
of the symplectic Chern number 


1. Berry phase of a Bogoliubov quasi-particle 


For an excitation conserving Hamiltonian, the Chern 
number of the ra-th band can be viewed as a sum of Berry 
phases accumulated on a set of closed loops covering the 
whole Brillouin zone. In this case, the state which accu¬ 
mulates the relevant Berry phase is the m-th eigenstate of 
the single-particle Hamiltonian h k (a block with quasimo¬ 
mentum k of the single-particle Hamiltonian h). Below, 
we show that one can naturally extend this definition of 
the Chern number to any bosonic Hamiltonian including 
anomalous terms by identifying the relevant Berry phase 
in a second-quantized setting. 

For each quasi-momentum k, the second-quantized 
block 


#k = 


a_ k 


M k) 


ak 

t 


a 


-k 


of the full bosonic Hamiltonian H is a six-mode squeez¬ 
ing Hamiltonian. If we regard the quasi-momentum k 
as an external parameter, we can ask ourself what is the 
additional Berry phase accumulated by a single Bogoli¬ 
ubov quasi-particle in a specific band n while the quasi¬ 
momentum is varied adiabatically over a closed loop. In 
other words, we calculate the Berry phase accumulated 
by the many-body state /?£ n |*Sk) where |5k) is the Bo¬ 
goliubov vacuum. We find 


Vn = *^<S k |/VVk/3£,JS k > • dk 

= iS(S k \\$^ n ,^ k ^JS k )-dk + id(S k \^ k \S k )-dk 

= i S (“k,nW V kMk,n[«] “ H Vk^k.n [«]) ' dk 

^(S k |V k |S k ) • dk 


= <j) A n - dk + i <^>(<S k |V k |S' k ) • dk. 


In the second line we have used that /3k,n|*Sk) = 0 (by 
definition of the vacuum). We note that the Bogoli¬ 
ubov vacuum 15k) is quasi-momentum dependend and 


could possibly accumulate a Berry phase by its own, 
i </>(Sk|Vk|Sk) • dk ^ 0. However, the Berry phase of in¬ 
terest is the additional Berry phase accumulated by the 
quasi-particle added over the Bogoliubov vacuum. 


2. Properties of the symplectic Chern numbers 

Taking into account the orthonormality condition 
Eq. one can immediately prove that the Chern num¬ 
bers have the usual properties: (i) They are integer num¬ 
bers; (ii) After a phase transition where two or more 
bands touch the individual Chern number of the band 
involved in the crossings may change but their sum does 
not change. Since the crossing of a particle and hole band 
lead to an instability rather than a phase transition, the 
sum of the Chern numbers over the particle bands is zero. 


Appendix C: Details of the calculation of the 
topological phase diagrams 

L. Symmetry of the topological phase diagram 
under synthetic magnetic field inversion 


In the topological phase diagram of Fig.[3^b) all Chern 
numbers change sign if the direction of the synthetic 
gauge field is inverted, —>• — 4>. This has a simple 
explanation: To change the sign of the flux <I> and of 
the quasimomentum k corresponds to taking the com¬ 
plex conjugate of the BdG Hamiltonian in momentum 

space, h(— k, — <f>) = h*( k, <f>), c.f. Eq. A3 It follows 
that the single-particle eigenfunctions for opposite val¬ 
ues of the flux and of the quasi-momentum are also re¬ 
lated by complex conjugation, |k n (4>)) = (| — k n (—<£>)))*. 
From the definition of the Chern numbers, c. f. Eq. [7j it 
immediately follows that the Chern numbers change sign 
under inversion of the synthetic gauge field <f>. 


2. Border of the different topological phases 

At a border of a topological phase transition a pair 
of Chern numbers can change their values because the 
corresponding bands touch. Generally speaking bands 
tend to repel each other rather than crossing. How¬ 
ever, at a lattice symmetry point this phenomenon does 
not necessarily occur because the interaction of a pair 
of bands can be prevented by a selection rule. In par¬ 
ticular, at the rotational symmetry points A, K\ and 
T, a hole with quasi-angular momentum m can only be 
converted into a particle with quasi-angular momentum 
m + m v . We note that due to inversion symmetry the 
bands must touch simultaneously at the symmetry points 
K and K '. We refer to the set of parameters where the 
bands touch at the symmetry points K and K' (T) as 
A-lines (T-lines). When also time-reversal symmetry is 
present there is band crossing at all rotational symmetry 
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points. We refer to the set of parameters where time- 
reversal symmetry occurs as T-lines. In addition, a pair 
of bands can touch at one of the three M points where 
one sublattice is decoupled from the remaining two sub¬ 
lattices (a particle or hole on that sublattice can not hop 
on the remaining sublattices). We note that due to ro¬ 
tational symmetry a pair of bands should touch simul¬ 
taneously at all three M points. We refer to the set of 
parameters where a pair of bands touch at the M points 
as M-lines. 

In our highly symmetric system, we expect most of the 
crossings to occur at a symmetry point. However, we note 
that accidental crossings away from any symmetry point 
are not forbidden. Indeed most (but not all) borders of 
the different topological phases in Fig. 2 can be identified 
with T-lines, T-lines, T-lines, or M-lines as explained 
below. 

We first focus on the effective model. The vertical lines 
<f> = 0 , ± 37 r are T-lines (there is time-reversal symmetry 
because the hopping amplitude J is real). One can also 
easily recognize the M-lines because they are horizontal. 
This must be the case because the spectrum at a M- 
point where a sublattice decouples from the remaining 
sublattices does not depend on the flux <f>. Indeed, there 
is such a horizontal line in the phase diagram of Fig. 2(b). 
We note that it appears for z > 2 ~ \Jij\&- Below, we show 
that this analytical expression holds when the a quasi¬ 
particles are described by an effective particle-conserving 
Hamiltonian. We can also find an analytical expression 
for the T-line and the T-lines as explained below. 

We initially focus on the T-lines. We regard the band 
crossing condition of a pair of levels with quasi-angular 
momentum m and m', T m (T) = T m /(T) as an implicit 
equation for the parametric coupling v as a function of 
the flux <f>. We take advantage of the analytical expres¬ 
sion for t he spectrum at the rotational symmetry points 
Eq. (A5) to solve this equation. For m = 0 and m' = — 1, 


we find exactly one real positive solution 


1 


v = 


2^3 


<f> 


UJq — 4 J cos — ) 2 — 


cjo — 4 J cos 


2tt ± < f > 


1 1/2 


in the intervals — 37 r < <f> < —it and 2tt < <f> < 37 r. 
Keeping in mind that the phase diagram is periodic with 
period 67 r, this solution can be thought of as a single T- 
line which goes (for increasing flux) from v = 0 at <F = 2ir 
back to v = 0 at <f> = 5tt (<T = — 7 r). Indeed, such a line 
is visible in the phase diagram of Fig. 2(b). From the 
implicit equations E m ( k) = E m /( k), k = T,T one can 
find similar formulas for the remaining T-line and the T- 
lines. In particular, the other T-line corresponds to the 
crossings of the levels with angular-momentum mm 1 
and m = 0 and goes from v = 0 at <f> = it back to v = 0 
at <F = 47 r (4> = —27r), see also Fig. 2(b). There is 
not a third T-line because the levels with quasi-angular 
momentum m — 1 and m — — 1 are degenerate only on 


the T-lines. Likewise, one can show that the T-lines go 
from v = 0 at 4> = —it back to v = 0 at <f> = it and from 
v = 0 at <f> = — 2tt back to v = 0 at <f> = 27 t, respectively. 
We note that the formulas for the band crossings are 
exact and valid for an arbitrary ratio of J/Co. However, 
if J/Co is above a finite threshold the unstable region may 
overlap with the band crossings and not all topological 
phases will be present in the phase diagram. 

Above we have identified all lines forming the border 
of the different topological phases in Fig. 2(b) except for 
the lines which appear above the M-lines very close to 
the T-lines and surrounds the white areas of the topolog¬ 
ical phase diagram Fig. 2 b. These lines correspond to 
accidental crossings which occur away from any symme¬ 
try point. They enclose four different topological phases 
(which are not listed in our legend for brevity). 

Next, we discuss the topological phase diagram of the 
original model. For a pump circulation m v = 1 (m v = 
— 1), the resulting effective flux <f> is positive (negative), 
see Eq. < 0 - Thus, topological phase diagram of the 
original model is a deformed version of the right (left) 
half of the effective model phase diagram, see panel (a) 
of Fig. 2 for the case m v = 1 (m u = —1). The case 
m v = 0 on the other hand is mapped onto the T-lines of 
the effective diagram. 

A remarkable feature of our model is that there is 
only a single topological phase for any fixed value m v 
of the pump circulation if the off-diagonal parametric 
terms are not present (i/ 0 q = 0): C = (=Fl,0, ±1) for 


m v — ±1 and C = ( 0 , 0 , 0 ) for m v = 0 . This is reminis¬ 
cent of the anomalous Quantum Hall effect on a Kagome 
lattice with nearest neighbor hoppings (OMN model) 
where the topological phase is uniquely determined by 
the sign of the magnetic flux piercing a triangular pla- 
quette, C = (=Fl, 0, ±1) if the flux is positive or negative, 
respectively. Indeed, for small squeezing u on , J <C cjo? 
the parametric interaction effectively induces a small syn¬ 
thetic gauge field with a positive flux for m u = 1. This 
can be easily seen by switching to the effective descrip¬ 
tion and neglecting the residual parametric terms. For 
concreteness we consider the case m v = 1. For small 
squeezing and = 0 we are somewhere close to v = 0 , 
<F = 2tt inside the topological phase C = (— 1 , 0 , 1 ) at 
the bottom right corner of the effective diagram. From 
the above analysis of the effective diagram we know that 
a topological phase transition can occur only if we cross 
a r or T-line. However, from the analytical solution of 
the spectrum at th e ro tational symmetry points of the 
original model Eq. (A4) we see that such crossings never 
occur on the v q r = 0 axis. Thus, there in no topological 
phase transition even for large squeezing if v q r = 0 . 


3. Effective excitation-conserving Hamiltonian 

When Co is much larger than \J\ and v one can derive 
an effective excitation conserving Hamiltonian. In the 
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regime where v \J\ it is not enough to keep the ex¬ 
citation conserving terms in Eq. (|8| but one should also 
include the leading order correction in v/ uj. We arrive at 
the excitation conserving Hamiltonian with next-nearest- 
neighbor hoppings, 


= E 4 2)& ] &i - ( C1 ) 

j <j,i> «j,i» 


Here, ((j,l)) indicates the sum over next-nearest- 
neighbor sites and 


2z> 




UJ = UJ — 


UJ 


7 ( 2 ) _ Yl 
^ _ 2uj 


j (!) _ f , ^ 


In this simplified picture it is straightforward to calcu¬ 
late the band structure at the M points and finding the 
bad degeneracy condition z> 2 = \Jj\\uJ which leads to the 
horizontal line in the topological phase diagram of the 
effective model. 


Appendix D: Bulk-boundary correspondence 


It is well known that in a system with a boundary, the 
net number of edge states (the number of right-movers 
minus the number of left-movers) in a bulk band gap 
is a topological invariant ES3- This statement is based 
on the sole assumption that the band structure and the 
corresponding eigenvectors change smoothly in the pres¬ 
ence of a local perturbation that does not close a gap. 
Thus, it clearly applies to any quadratic Hamiltonian. 
For the special case of an excitation-conserving insulator 
or a superconductor the bulk boundary correspondence 
expresses such topological invariant in terms of the Chern 
numbers: the net number of edge states in a band gap 
coincides with the sum of the Chern numbers of all bands 
below that band gap. Here, we explicitly show that the 
bulk-boundary correspondence is still valid for our model 
where anomalous terms are present. 

We start noticin g th at the wavefunctions of the RWA 
Hamiltonian Eq. (Cl) depend only on two parameters: 


the phase <f>/3 of J and the dimensionless next-nearest- 
neighbor coupling z> 2 /(| J\uj). By calculating the phase di¬ 
agram as a function of these parameters [not shown] we 
see that it supports all topological phases present in the 
topological phase diagram of the effective Hamiltonian 
for the d-quasiparticles [our full model without approxi¬ 
mations]. Thus, we can continously interpolate between 
the two Hamiltonians without crossing any topological 
phase transition [by sending uj —> oc while also tuning 
<f> and i> 2 /(uj\J\) to stay in the same topological phase]. 
Keeping in mind that the bulk-boundary correspond ence 
holds for the excitation-conserving Hamiltonian Eq. ( |C1| ) 
and that the net number of edge states does not change 
during the interpolation [unless a gap is closed], we can 


conclude that such correspondence is valid for our mo del 
even for small uj where the RWA leading to Eq. (Cl) is 
not a good approximation. 

We note that the above reasoning combined with the 
assumption that a continuous interpolation between any 
quadratic bosonic Hamiltonian and an excitation con¬ 
serving Hamiltonian is always possible without closing 
any band gap, leads to the general validity of the bulk- 
edge correspondence. 


Appendix E: Details of the transport calculations 

In our transport calculations we have included photon 
decay. We adopt the standard description of the dis¬ 
sipative dynamics of photonic systems in terms of the 
Langevin equation of input-output theory |33], for each 
site: 

aj = dj\ — Kctj /2 + (El) 

In practice, we consider an array of detuned parametric 
amplifiers with intensity decay rate k and add to the stan¬ 
dard description of each parametric amplifier the inter¬ 
cell coherent coupling described in the main text. The 
last term describes the input field and includes the vac¬ 
uum fluctuations as well as the influence of an additional 
probe field. The field leaking out of each cavity at site j is 
given by the input-output relations a)° ut ^ = a^ — y/ddj . 
The above formulas give an accurate description of a pho¬ 
tonic system where the intrinsic losses during injection 
and inside the system are negligible. 

In Fig. (3), we show the probabilities Te(cj,/,j) and 
T/(cj, /, j) that a photon injected on site j with frequency 
cjin = cj + co’l/ 2 is transmitted to site l elastically (at 
frequency uj+ujl/2) or inelastically (at frequency ujl/2 — 
uj) where it is detected. From the Kubo formula and the 
input-output relations we find 

Te(W: l,j) = |Sij - l,j) | 2 

T I (LU,l,j)=K, 2 \G I (uj,l,j)\ 2 . 

They depend on the Green’s function in frequency space 
G(u,l,j) = dte iut G(t, l,j) where 

,._(G E (t,l,j) Gj(t, l,j) 
hJ >~ \G T (t,l,j) G* E (t,l,j) 

with the elastic and inelastic components 

G E {t,l,j) = -iQ(t)([ai(t),a t(0)]) and Gi(t,l,j) = 
— i@(i)([aj(i),aj(0)]), respectively. From the Kubo for¬ 
mula and the input-output relations a.- '" 1 ' = a - n) — Jncij 
we find 

T E (u>,l,j) = |% -mG E (u,l,j)\ 2 (E2) 

T I (cu,l,j) = K 2 \G I (u;,l,j)\ 2 . (E3) 

Taking into account that 

&j(t) = E e- iE ^u n [j}b n + e iE ^v* n \jp n 
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where b n are the ladder operators of the normal 
modes for a finite array of N sites and | n) = 
( u n [1],..., u n [TV] , v n [1],..., v n [TV] ) T are the correspond¬ 
ing single-particle states, the Green’s function reads 


Ge(u, l,j) 


u n\l\Ki\j\ 

^ uo — E[n\ + in /2 


Gi(u,l,j) 


v n [%n[i] 

^ uo — E[n] + in /2 


uo -|- £/[??/] in/‘l 
(E4) 

u ni l ] v n\j] 

uo J- -£/[?7<] T in/ 2 

(E5) 


We note that for a probe field inside the bandwidth of the 
particle (hole) sector but far detuned from the hole (par¬ 
ticle) sec tor, only t he fi rst (second) term of the summand 
in Eq. (|E4|) and (E5) is resonant. Thus, as expected, 


the inelastic scattering is comparatively larger when the 
probe field is in the hole band gap. 

It is easy to estimate quantitatively the relative in¬ 


tensities of elastically and inelastically transmitted light 
when the parametric interaction of the a quasiparti¬ 
cles is small. In this case, it is straighforward to show 
that | v n [j ] / [j ] | ~ tanhr independent of the band n 
and the site j. By putting together Eqs. flE2|E3|E4|E5| 
and neglecting the off-resonant terms we find that for 
uo,uo > | J|, n : | uo — uo |, 

Ti(u),l,j) « (tanh r) 2 T E (uJ,l,j) 

« Ti(-u,l,j) « (coth r) 2 T E (-uj,l,j). 


These analytical formulas agree quantitatively with the 
numerical results shown in Fig. 4(b) [note that in Fig. 
4(b) the transmission at the output sites is rescaled by 
the overall transmission, Ti(u,l, j) + Te(uo, /, j)]. 

Note that we have assumed, for simplicity, that there is 
no intrinsic absorption present. If there is intrinsic pho¬ 
ton absorption, that will add another decay channel to 
the equation for the light field, but the resulting picture 
for the light field propagation remains unchanged except 
for the expected reduction in propagation length along 
the edge state. 
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